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We study inelastic cotunneling through a strong Coulomb-blockaded quantum dot subject to 
a static magnetic field and a perpendicular circularly-polarized magnetic field using a quantum 
Langevin equation approach. Our calculation predicts an interesting controllable cotunneling cur- 
rent characteristic, splitting-zero-anomaly-splitting transition of the differential conductance with 
increasing the driving frequency, ascribing to the role of photon- assisted spin-flip cotunneling pro- 
cesses. 
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Coherent control of spin dynamics in a quantum dot 
(QD) using circularly-polarized magnetic field (CPF)^'^ 
in a solid-state environment has been a subject of in- 
creasing interest in recent years 4 This interest has led 
to a large amount of work on detection of the electron 
spin resonance (ESR) in a QD)^!^!^!^ and initializing the 
state of electron spin via optical absorption of circularly 
polarized light®'^''^°. Recently, the single-electron ESR in 
a QD-lead system with sizable Zeeman splitting has been 
theoretically reported to carry pure spin flow, which can 
be used as a fundamental element and/or a spin source 
device in the spintronic circuit. ^^'^^ These papers have 
studied the transport through a QD in resonant tunneling 
regime by the nonequilibrium Green's function (NGF)^ 
and the quantum rate equation approachfiS, respectively. 
The effect of the s-d exchange interaction between the lo- 
calized spin and conduction electrons on electron trans- 
mission probability in linear transport regime has been 
exploited for a mesoscopic ring embedded with an mag- 
netic impurity by the NGF.^'^ In this paper, we focus 
our studies, for the first time, on the influence of a driv- 
ing CPF upon cotunneling (non-resonant coherent tun- 
neling) through a strongly interacting QD in the weak- 
tunneling limit fi^ii^ predicting that the cotunneling cur- 
rent can be easily controlled by tuning the driving fre- 
quency. 

Cotunneling through a single-level Cd interacting QD 
subject to an ambient constant magnetic field Bq along 
the z-axis and a driving external (magnetic) CPF with 
frequency Uc whose direction rotates in the plane perpen- 
dicular to the z-axis, Bi(i) — i3i(cosLi;ci,sinLi;ci, 0), can 
be described by the Hamiltonian H = Hb + Hq + Hi'}'^ 
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where 9 — ujct, cj^ko- (^Tykcr) is the creation (annihilation) 
operator for electrons with momentum k, spin-cr and en- 
ergy e^k in lead (= L,R), S = (S*^, 5'?', 5^) are Pauh 
spin operators of electrons in the QD [S^ = S*^ ± iS^], 
and Jrirji is the exchange coupling constant. Aq = q^bBq 
is the static magnetic-field So-induced Zeeman term, and 
A = g^BBi describes the spin-flip scattering caused by 
the CPF. 7?dir is the potential scattering term, which is 
decoupled from the electron spin due to the number of 
electrons in the dot level being one. As a result, this 
term has no influence on the dynamical evolution of the 
electron spin and behaves only as a direct bridge to con- 
nect the left and right leads. For an Anderson model 
with symmetrical coupling to the leads t, we have J^ri' — 
2Jo = I td- As in our previous paper—, we can rewrite 
the tunneling term, Eq. \\c\ . as a sum of three products 
of two variables: Hi = Q^S*^ -h Q+S'" + + with 

the same definitions of Q^*-*-* as in Ref. 16 and = Hdh-- 
It is noted that our model describes the non-resonant 
tunneling through a localized magnetic impurity involv- 
ing s-d exchange interaction with the conduction elec- 
trons in the presence of both a static magnetic field Bp 
and a rotating magnetic field Bi. We assume the ex- 
change interaction to be weak so that no Kondo effect 
emerges, and assume that charge fluctuation completely 
vanishes. 

It is well-known that the isolated spin- 1/2 electron 
under the influence of both Bq and an external CPF 
Bi(i) (the Rabi problen*^), Eq. (jlbp . has an analytical 



solution. ' In the rotating frame, x — cos 6*5^ — sin 



y = sin6l5^ +cos6l5'^ and z 
Eq. (jlbp takes the form 



, the QD Hamiltonian 



Ho 



-Ax — 5z, 
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with 6 — Ao — ujc- Correspondingly, the Heisenberg equa- 
tions of motion (EOM's) of these free rotating coordi- 
nates become: x = Sy, y = —Sx -f Az, and i = —Ay. 
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Solving these resulting ordinary differential equations, we 
obtain the free evolutions of these rotating coordinates: 



r{t') = M(r)r(i), [r{t) ^ {x{t),y{t), z{t))' 
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frequency fl = VA^ + S'^. 

Furthermore, the transformed interacting Hamilto- 
nian, Eq. (jlcp . reads in terms of these rotating frame: 

Hi = Q''x + Qyy + Q^z + Q\ (4) 

with = e-'''g~+e*^Q+ and Q?' = i{e-'^Q- -e''^Q+). 
The Heisenberg EOM's for the spin operators are: 



y = -5x + Az + Q^'x - Q^'z, 
i = -Ay-Qyx + Q^'y. 



(5a) 
(5b) 
(5c) 



It is clear that the spin dynamics, apart from free evo- 
lutions, are perturbatively modified by the weak tunnel 
coupling. To obtain the modified dynamics, we employ a 
generic quantum Langevin equation approach. In 
our derivation, operators of the QD spin and the reser- 
voirs are first expressed formally by integration of their 
Heisenberg EOM's, Eq. (O, exactly to all orders of J^^'. 
Next, under the assumption that the time scale of decay 
processes is much slower than that of free evolutions, we 
replace the time-dependent operators involved in the in- 
tegrals of these EOM's approximately in terms of their 
free evolutions, Eq. Thirdly, these EOM's are ex- 
panded in powers of Jjjrj' up to second order. To this end, 
we can establish the Bloch-type dynamical equations for 
the averaged spin variables r as: 
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with the reservoir correlation function, C{lj)^ and the 
response function, R{uj), defined as 

{uj + V\ foj-V'^' 
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with griri' — J^rj'Po (Po ^^le constant density of states 
of both electrodes) and (p{x) = a;coth(x/2). V is the 
bias- voltage applied symmetrically to the two electrodes, 
/iL = —PR — V/2. We use units with h — ks = e = 1. 

Notice that these quantities of Tap describe the dissi- 
pation of cotunneling processes to the dynamics of the 
QD spin, in conjunction with the effect of CPF: transi- 
tion between two spin-splitting states ±17 /2 (spin analog 
of charge optical Stark effect) due to photo-absorption 
and/or emission. Also, it is easy to see that the Bloch 
equations, Eq. ([5]), is exactly reduced to the results of 
Ref. [l^ in the case of vanishing CPF, A = and Wc = 0. 

The nonequilibrium steady-state spin projections of 
the QD in the rotating frame, y°° = {x°° ,y°° , z°°), 
can be readily obtained from Eq. ©. As is evident, 
x°° = y°° = 0, while z°° = ii?(Ao)/C(Ao) is identi- 
cal to the previous theoretical result of nonequilibrium 
magnetizationi^iSS in absence of driving CPF. On the 
other hand, in the case of vanishing static magnetic field, 
Aq = 0, the nonzero driving CPF can induce an ad- 
ditional spin orientation z°° ^ along the rotating di- 
rection of the CPF, together with the nonzero 2:°° (not 
shown here). At equilibrium, this optically-induced spin 
orientation phenomenon has historically been termed, in 
literature, as either the inverse Faraday effect^ or as 
the Zeeman light shift, which is ascribed to the CPF- 
induced spin-splitting (ac spin Stark effect). In addi- 
(7c) tion, one can interestingly observe a nearly vanishing y- 
component of the spin polarization even at CPF driving. 

For nonzero static magnetic field, previous theoretical 
studies show that coherent suppression of tunneling may 




take place only if the Rabi frequency Q matches with the 
driving frequency ujc, which corresponds to the following 
relationship of the driving field^i: 



UJr 



2Ao 
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In Fig. 1, we plot the nonequilibrium spin projection 
under the influence of both nonzero static magnetic 
field and driving CPFs. Different from the results of 
Aq = 0, we find an obvious nonzero polarization of the y- 
component of QD spin near the resonant field, Wc = Aq, 
for the weak driving field A ~ 0.2Ao [Fig. 1(a)]. Inter- 
estingly, we observe z°° = at = Aq for any strength 
of CPF, meaning that the resonant frequency of a CPF 
may overcome the Zeeman-splitting of a static magnetic 
field. At the same time, x°° becomes unpolarized for 
weak driving field [Fig. 1(a)], while reaches the maxi- 
mum value for strong driving field, which decreases with 
increasing of the bias- voltage [Fig. 1(b)]. Moreover, we 
notice that (1) a:°° = = at = 20*; and (2) z°° is 
negative if Ag < Wc < 2r2* for strong driving CPF. 

We now proceed with the calculation of tunneling 
current. The current operator through the QD is de- 
fined as the time rate of change of charge density = 

a '^?)k(T'^»)k(T 1^9''^ 'n ('^6 choose the left lead as an 
example) : 



+^iQlR - Qrl)s+ + ^{QIr - QtL)S~ 
- Q^l + Q^R - (10) 



where the definitions of Q^^r and Q^^, can be found in 
our previous paper— and those terms in the last line 
are stemming from the direct tunneling term i?dir — 
Jo J2riri'.a Qfj'^'' ■ Thc fincar-respouse theory gives 



i = {Mt)) 



dt'{[JL{t),Hi{t')Ua, (11) 



where the statistical average (• ■ ■ )o is performed with re- 
spect to decoupled two subsystems, QD and reservoirs. 
Inserting Eqs. (fTc|) and PH)) into Eq. (fTT|) . one can derive 
the explicit expression for steady-state current in terms 
of the steady-state spin projections r°° through a lengthy 
but straightforward calculation. In particular, we note 
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because the QD is connected to two normal leads in the 
system under consideration^^ These results indicate that 
the contribution of the direct tunneling term to the cur- 
rent is independent of the dynamics of the QD and only 
depends on the bias-voltage and temperature of the two 
electrodes. To this end, the current / is 
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Different from the previous resuhs of the nondriving 
QD system/® the nonzero a;-component of the stationary 
spin polarization has additional contribution to the cur- 
rent besides the z-component. The linearly bias- voltage- 
dependent term in Eq. p6p stems from the cotunnel- 
ing processes, in which the spin projection remain un- 
changed. The CPF has no effect on these processes. 
More importantly, one observes that the current formula 
Eq. (|16p includes some new factors involving uJc±U,±V 
[Eq. P?|) ]. which can be intuitively ascribed to the photo- 
assisted spin-flip cotunneling processes. For instance, the 
term involving — cjc results from the contribution of co- 
tunneling accompanied by an absorption of one-photon 
[Fig. 2(a)], while the term involving + lOc corresponds 
to the process assisted by a spontaneous emission of one- 
photon [Fig. 2(b)]. Therefore, we simply name these co- 
tunnelings as SI ± processes respectively below. Be- 
sides, the factors Q,±V m. Eq. stem from the spin- 
flip cotunneling events without spontaneous emission or 
absorption of photon, while the role of driving field is 
reflected through the effective spin-splitting SI. In the 
following, we will see that the joint effects of these new 
terms are responsible for the controllable patterns of co- 
tunneling through a QD by the presence of CPF. 
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FIG. 2: Schematic diagrams of photon-absorption- {Q, — loc) 
(a) and photon-emission- {Q, + ujc) (b) assisted spin-flip co- 
tunneling processes. 



We plot the differential conductance G — ^ (in unit 
of 2ngLji) vs bias- voltage as functions of different driv- 
ing frequencies ujc for a weak driving CPF A/Aq — 0.2 



(a,b) and a strong driving field A/Aq = V2 (c,d) in 
Fig. [31 li ojc — 0, the I-V curve reduces to the ordinary 
cotunneling characteristics under an ambient constant 
magnetic field:^® a jump at = ±Q — zb-y/Ap + A2 
Ao in Fig. [3^a) due to A ^ Aq and — VSAq in 
Fig. [Sl^c), respectively], ascribing to the energetic inacti- 
vation (activation) of spin-flip processes at bias-voltage, 
\V\ < (>)S^- While for nonzero driving frequency, such 
as Wc/Aq = 0.5, the CPF effectively suppresses spin- 
splitting Sl/Ag « 0.54 at the case of weak driving held, 
thus leading to two jumps, one of which is located at 
V — ±r2 due to the pure activated spin-flip process, and 
another of which appears at V = zt{n + ujc) ~ ±Ao, 
which can be ascribed to opening of an additional channel 
for electron transfer cotunneling due to the one-photon- 
emission-assisted spin-flip resonant event. In contrast, 
the strong driving field A = \/2Ao makes S2 = 1.5Ao 
and generates three jumps: the first one corresponds to 
the Vl — ujc process; the second one is due to the pure elec- 
tronic spin-flip event (SI); and the third one results from 
the fi-l-Wc process. When the frequency increases to a bit 
higher than the resonant frequency [I.OAq (SI* = 1.5Ao) 
for the weak (strong) driving held], a zero-bias peak 
emerges for both cases, which can be understood by the 
fact that the proper high frequency CPF has enough en- 
ergy to spur the spin-flip cotunneling even at equilib- 
rium [see Fig. 2(a)], while the rising voltage, V > ±S7 
[Figs. M.a., b)] or ±(S7 - ujc) [Figs. M.C, d)], contrarily 
suppresses its activation until a pure bias-driven spin-flip 
process is excited a.tV — ±Ao or ±S1. Recently, the zero- 
anomaly (ZA) behavior of cotunneling current has been 
also reported for a QD connected to two anti-parallel fer- 
romagnetic electrodes. 25-26 Here, it should be pointed out 
that the appearance of ZA in the present system is due 
to photon-assisted inelastic spin-flip scattering becoming 
resonant in the presence of external magnetic field, which 
is different from the previous mechanism, elastic spin-flip 
event in the case of polarization leads without any mag- 
netic field ]25i2® As a result, the present system has more 
rich transport features by tuning strength and frequency 
of driving CPF. 

For the weak CPF field with uJc ~ 2Ao, the photon- 
assisted excitation of spin-flip cotunneling merges coinci- 
dently into pure bias excitation, leading to the disappear- 
ance of Z A behavior [Fig. [3{b)] . While for strong CPF 
field [Fig. [3fd)], the ZA also vanishes if > 2S7*, which 
is due to the peculiar features of the nonequilibrium spin 
projections and z°° [see Fig. 1(b)]. For higher driving 
frequency, the cotunneling exhibits splitting diflerential 
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FIG. 3: The calculated differential conductance dl/dV vs. bias-voltage V/Ao (Ao = 1.0) under weak driving CPF A/Ao = 0.2 

(a,b) and strong CPF A/Ao = \/2 (c,d) with several driven frequencies at a nonzero static magnetic field. The arrow indicates 
the direction of increase of driven frequency uic- The other parameters are the same as in Fig. 1. 



conductance and eventually tends to the pattern in the 
presence of a static magnetic field alone, because the QD 
spin cannot, from physical point of view, catch the de- 
tails of the driving field with considerably high frequency. 
Finally, we point out that the ZA feature is robust over 
a wide region of temperature (not shown here). 

Conclusion. — In summary, we have presented an ana- 
lytical study of the inelastic cotunneling, including the 
nonequilibrium spin projections and currents, in a sin- 



gle QD under a static magnetic field and a perpendicular 
CPF, revealing a controllable I-V pattern, the transition 
between ZA and splitting of the differential conductance, 
due to photon- assisted spin-flip inelastic cotunneling. 
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